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Abstract 

The process ^(3686) — > Tr°h c , h c — > ^r\ c has been studied with a data sample of 106 ± 4 million 
V>(3686) events collected with the BESIII detector at the BEPCII storage ring. The mass and 
width of the P-wave charmonium spin-singlet state h c (}P\) are determined by simultaneously 
fitting distributions of the 7r° recoil mass for 16 exclusive r\ c decay modes. The results, M(h c ) = 
3525.31±0.11 (stat.) ±0.14 (syst.) MeV/c 2 and F(h c ) = 0.70±0.28±0.22 MeV, are consistent with 
and more precise than previous measurements. We also determine the branching ratios for the 16 
exclusive r) c decay modes, five of which have not been measured previously. New measurements 
of the rj c line-shape parameters in the El transition h c — > jr/ c are made by selecting candidates 
in the h c signal sample and simultaneously fitting the hadronic mass spectra for the 16 rj c decay 
channels. The resulting r\ c mass and width values are M(rj c ) = 2984.49 ± 1.16 ± 0.52MeV/c 2 and 
r(r/ c ) = 36.4 ± 3.2 ± 1.7 MeV. 

PACS numbers: 14.40.Pq, 13.25. Gv, 12.38. Qk 
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I. INTRODUCTION 



Studies of charmonium states have played an important role in understanding Quantum 
Chromodynamics (QCD) because of their relative immunity from complications like rela- 
tivistic effects and the large value of the strong coupling constant a s . In the QCD potential 
model [H , the spin- independent one-gluon exchange part of the cc interaction has been de- 
fined quite well by existing experimental data. The spin dependence of the cc potential is 
not as well understood. Until recently, the only well-measured hyperfine splitting was that 
for the IS states of charmonium, AMhAlS) = MU/^) - M(t] c ) = 116 ± lMeV/c 2 Q. In 
the past several years BaBar Q, Belle 0, CLEO Q, and BESIII Q have succeeded in iden- 
tifying r) c (2S) and have measured AM h£ (2S) = M(^(3686)) - M(r) c (2S)) = 47 ±1 MeV/c 2 . 

Of the charmonium states below DD threshold, the /i c (l 1 Pi) is experimentally the least 
accessible because it cannot be produced directly in e + e~ annihilation or in the electric- 
dipole transition of a J PC = 1 charmonium state. Limited statistics and photon-detection 
challenges also were major obstacles to the observation of h c in charmonium transitions. The 
precise measurement of h c properties is important because a comparison of its mass with the 
masses of the 3P states (Xcj) provides much-needed information about the spin dependence 
of the cc interaction. According to QCD potential models, a small spin-spin interaction 
would lead to hyperfine splitting of the P states that is close to zero 0]. 

The first evidence of the h c state was reported by the Fermilab E760 experiment 
and was based on the process pp — > ir J/ip. This result was subsequently excluded by 
the successor experiment E835 [9], which investigated the same reaction with a larger data 
sample. E835 also studied pp — >■ h c — » 77/c in this case finding an h r signal. Soon after this 
the CLEO collaboration observed the h c and measured its mass ll[ by studying the 
decay chain 0(3686) — > n°h c ,h c — > yrj c in e + e~ collisions. CLEO subsequently presented 
evidence for h c decays to multi-pion final states [12]. Recently, the BESIII collaboration 
used inclusive methods to make the first measurements of the absolute branching ratios 
#0(3686) n°h c ) = (8.4 ± 1.3 ± 1.0) x 10" 4 and B(h c 7 t? c ) = (54.3 ± 6.7 ± 5.2)% 0. 
CLEO has confirmed the BESIII results |l4J] and also observed h c in e + e~ — > n + ir~h c at 



y/s = 4170 Me V, demonstrating a new prolific source of h c [15 . 

77 C (1S') is the lowest-lying ,5-wave spin-singlet charmonium state. Although it has been 
known for about thirty years [l6|, its resonance parameters are still interesting. For a 
long time, measurements of the r] c width from B-factories and from charmonium transitions 
were inconsistent . The discrepancies can be attributed to poor statistics and inadequate 
consideration of interference between t] c decays and non-resonant backgrounds. Besides, 
the rj c line shape also could be distorted by photon energy dependence in the Ml(or El) 
transition, which will affect the resonance parameters measurement. Recent studies by 



Belle, Babar, CLEO, and BESIII 17H20|. with large data samples and careful consideration 
of interference, obtained similar 7] c width and mass results in two-photon-fusion production 
and -0(3686) decays. The h c — > ^rj c transition can provide a new laboratory to study r\ c 
properties. The t] c line shape in the El transition h c — > 7?7 C should not be as distorted as in 
other charmonium decays, because non-resonant interfering backgrounds to the dominant 
transition are small. 

In this paper, we report new measurements of the mass and width of the h c and i] c , and of 
the branching ratios i?i(^(3686) — > n°h c ) x i3 2 (^ c — > IVc) x ^(Vc Aj) and B 3 (r] c — > Aj), 
via the sequential process -0(3686) — > n°h c , h c — > •yr/c, rj c — > Aj. In this reaction Aj signifies 
16 exclusive hadronic final states: pp, 2(7r + 7r~), 2(A + A~), K + K~tt + ti~ , ppn + 'K^ , 3(7r + 7r~), 



4 



7r + 7r _ 7r°7r°, and 2(7r + 7r~)7T 7r°. Here Kg is reconstructed in its tt + tt~ decays, and 77 in its 77 
final state. The data sample of -0(3686) events was collected with the BESIII detector at 
the BEPCII e + e~ storage ring. 

The remainder of this paper is structured as follows: Sect. [TT] describes the experiment 
and data sample; Sect. [Ill] presents the event selection and background analysis; Sect. HVl 
discusses the extraction of h c and r\ c results; Sect. [V] describes the estimation of systematic 
uncertainties; and Sec. [VI] provides a summary and discussion of the results. 



II. EXPERIMENT AND DATA SAMPLE 



BEPCII is a two-ring e + e~ collider designed for a peak luminosity of 10 33 cm~ 2 s _1 at a 
beam current of 0.93 A per beam. The cylindrical core of the BESIII detector consists of 
a helium-gas-based drift chamber (MDC) for charged-particle tracking and particle identi- 
fication by dE/dx, a plastic scintillator time-of-flight system (TOF) for additional particle 
identification, and a 6240-crystal CsI(Tl) Electromagnetic Calorimeter (EMC) for electron 
identification and photon detection. These components are all enclosed in a superconduct- 
ing solenoidal magnet providing a 1.0- T magnetic field. The solenoid is supported by an 
octagonal flux-return yoke with resistive-plate-counter muon detector modules (MUC) in- 
terleaved with steel. The geometrical acceptance for charged tracks and photons is 93% of 
Att, and the resolutions for charged-track momentum and photon energy at 1 GeV are 0.5% 
and 2.5%, respectively. More details on the features and capabilities of BESIII are provided 
in Ref. jS). 

The data sample for this analysis consists of 156.4 pb _1 of e + e~ annihilation data collected 
at a center-of-mass energy of 3.686 GeV, the peak of the -0(3686) resonance. By measuring 
the production of multihadronic events we determine the number of 0(3686) decays in the 
sample to be (1.06 ±0.04) x 10 8 , where the uncertainty is dominated by systematics [22]. An 
additional 42 pb _1 of data were collected at a center-of-mass energy of 3.65 GeV to determine 
non-resonant continuum background contributions. 

The optimization of the event selection and the estimation of physics backgrounds are 
performed with simulated Monte Carlo (MC) samples. A GEANT4-based (231 . 24 1 detector sim- 
ulation package is used to model the detector response. Signal and background processes are 



generated with specialized models that have been packaged and customized for BESIII [25 



The 0(3686) resonance is generated by KKMC [26], and EvtGen [27( is used to model events 
for -0(3686) — > TT°h c and for exclusive backgrounds in -0(3686) decays. An inclusive sample 
(100 million events) is used to simulate hadronic background processes. Known -0(3686) 
decay modes are generated with EvtGen, using branching ratios set to world-average val- 
ues Q. The remaining 0(3686) decay modes are generated by LUNDCHARM j25|, which is 
based on JETSET [28] and tuned for the charm-energy region. The decays -0(3686) — > ir°h c 
are excluded from this sample. 

The 0(3686) — > n°h c events are generated with an h c mass of 3525.28 MeV/c 2 and a 
width equal to that of the xa (0.9 Me V). The El transition h c — > / -fi] c is generated with 
an angular distribution in the h c rest frame of 1 + cos 2 #*, where 9* is the angle of the El 
photon with respect to the beam direction in the h c rest frame. Multi-body r] c decays are 
generated according to phase space. 
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III. EVENT SELECTION AND BACKGROUND ANALYSIS 



For ^(3686) — > ir h c , h c — > / ~fr] c , the expected 7r° momentum is P n o ~ 84MeV/c, and the 
El transition photon emitted in h c — > jrjc has an expected energy of E^ei) — 503 MeV 
in the h c rest frame. We select signal candidates by demanding consistency with these 
expectations. For the selected candidates we then compute the n° recoil mass. Fitting the 
distribution of this recoil mass for the full event sample gives the results for the h c resonance 
parameters and signal yields. 

Charged tracks in BESIII are reconstructed from MDC hits within a polar-angle (9) 
acceptance range of |cos#| < 0.93. To optimize the momentum measurement, we require 
that these tracks be reconstructed to pass within 10 cm of the interaction point in the 
beam direction and within 1 cm in the plane perpendicular to the beam. Tracks used in 
reconstructing Kg decays are exempted from these requirements. 

A vertex fit constrains charged tracks to a common production vertex, which is updated 
on a run-by-run basis. For each charged track, time-of-flight and dE/dx information is 
combined to compute particle identification (PID) confidence levels for the pion, kaon, and 
proton hypotheses. The track is assigned to the particle type with the highest confidence 
level. 

Electromagnetic showers are reconstructed by clustering EMC crystal energies. Efficiency 
and energy resolution are improved by including energy deposits in nearby TOF counters. 
A photon candidate is defined as a shower with an energy deposit of at least 25 MeV in 
the "barrel" region (| cos 6*1 < 0.8), or of at least 50 MeV in the "end-cap" region (0.86 < 
| cos 6*1 < 0.92). Showers at angles intermediate between the barrel and the end-cap are 
not well measured and are rejected. An additional requirement on the EMC hit timing 
suppresses electronic noise and energy deposits unrelated to the event. 

A candidate 7r°(?7) is reconstructed from pairs of photons with an invariant mass in the 
range |M 77 - m n o\ < 15MeV/c 2 (|M 77 - m v \ < 15MeV/c 2 ) @. A one-constraint (1-C) 
kinematic fit is performed to improve the energy resolution, with the M{ r y r y) constrained to 
the known 7r°(?7) mass. 

We reconstruct K s — > 7r + 7r~ candidates using pairs of oppositely charged tracks with an 
invariant mass in the range \M nw — m K o | < 20 MeV/c 2 , where m K o is the known Kg mass 0]. 
To reject random tt + tt~ combinations, a secondary-vertex fitting algorithm is employed to 
impose the kinematic constraint between the production and decay vertices j29[. Accepted 
K s candidates are required to have a decay length of at least twice the vertex resolution. 

The i] c candidate is reconstructed in 16 exclusive decay modes, and the event is ac- 
cepted or rejected based on consistency with the h c — > •yrjc hypothesis. Specifically, the 
reconstructed mass M(rj c ) is required to be between 2.900 GeV/c 2 and 3.050 GeV/c 2 , and 
the transition-photon energy is required to be between 0.450 GeV and 0.550 GeV. Events 
passing this selection are subjected to a 4-constraint (4-C) kinematic fit to take advantage 
of energy- momentum conservation. Because of differing signal/background characteristics, 
we individually optimize requirements on x\c-> the X 2 °f the 4-C fit, for the 16 r] c chan- 
nels. If multiple T] c candidates are found in an event, the one with the smallest value of 
X 2 = Xlc + Xic + X 2 P id + Xvcrtcx is accepted, where xlc is the X 2 of the 1-C fit of the vr ^), 
Xp id is the PID \ 2 summation for all charged tracks included in the h c candidate, and Xvertex 
is the x 2 °f the K s vertex fit. If there is no 7r°/r] (Kg) in an event, the corresponding \\c 
(Xvertex) set to zero. 
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Based on studies of the inclusive MC sample, we identified several background processes 
with potential to reduce the precision of measurements made with specific r] c exclusive chan- 
nels because of sizable low-energy 7r° production. The processes and suppression procedures 
are as follows: 

• ^(3686) ->■ 7T+7T- J/il> 

The mass Mj of the system recoiling against the tt + tc~ in ^(3686) — > 7t + tt~X is 
calculated and the candidate is rejected if Mx is within ±12MeV/c 2 of the known 
J /if; mass. 

• ^(3686) -»■ n°7i°J/ip 

The mass Mx of the system recoiling against the 7r°7r° in ^(3686) — > 7r°7r°X is cal- 
culated and the candidate is rejected if Mx is within ±15MeV/c 2 of the known J /if} 
mass for all f] c final states except 7r + 7r~7r + 7r~7r°7r°. For this mode the lower tt° mo- 
mentum leads to recoil masses near 3.1GeV/c 2 , so the exclusion window is narrowed 
to±10MeV/c 2 . 

• ^(3686) -> 7x c2 

A candidate is rejected if it includes a 7r° for which either daughter photon has an 
energy within ±5MeV of that expected for the ^(3686) radiative transition to Xc2 
(128 MeV). 

• El photon candidates that are tt° decay products 

A candidate is rejected if its El photon can be combined with another photon in the 
event to form a 7r° within a mass window of ±10MeV/c 2 . 

• 7T° candidates that are from r) — > 7r + 7r~7r° 

Masses M(7r + 7r~7r°) are calculated for all possible combinations in the event and the 
candidate is rejected if any combination has a mass within ±15MeV/c 2 of the known 
rj mass. 

Decisions about whether to apply a requirement to a particular r\ c mode and the opti- 
mization of the x!c an d PID requirements were made on a channel-by-channel basis. The 
figure-of- merit used was S = N$ / \J N$ + Nb, where N$ is the number of signal and iVg 
the number of background candidates. PDG values [2| are used for the input rj c branching 
ratios, and for channels not tabulated by the PDG we estimate branching ratios based on 
conjugate channels or other similar modes. The optimized selection criteria are listed in 
Table [J in which the N(p), N(tt) and N(K) denote the numbers of identified protons, pions 
and kaons in a event. 

The 7T° recoil mass spectra for events passing these requirements show clear h c signals in 
the expected range, as can be seen in Fig. [U No peaking backgrounds in the signal region 
are found in the 100-million-event inclusive MC sample, in the continuum data sample taken 
at y/s = 3.65 GeV, or in ?7 c -candidate-mass side-band distributions. 
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TABLE I: Event-selection requirements for each exclusive channel. 
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FIG. 1: The vr° recoil mass spectrum in xp(3686) — > ir°h c , h c — > jrjc, rj c — > Xi summed over the 16 
final states Xj. The dots with error bars represent the 7r° recoil mass spectrum in data. The solid 
line shows the total fit function and the dashed line is the background component of the fit. 



IV. EXTRACTION OF YIELDS AND RESONANCE PARAMETERS 

We obtain the h c mass, width and branching ratios from simultaneous fits to the 7r° recoil 
mass distributions for the 16 exclusive r] c decay modes. Here only 1-C kinematic fits with 
7r° mass hypothesis are used to improve the energy resolution. The 4C-frts used in event 
selection are not used in the 7r° recoil mass reconstruction, because the energy resolution of 
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the signal ir° in 4C-fits is not as good as in the lC-fits, according to a MC study. From the 
same data sample we also determine the r\ c resonance parameters by fitting the 16 invariant- 
mass spectra of the hadronic system accompanying the transition photon in h c — > ^rjc. 

A. Fitting the h c signal 

To extract the h c resonant parameters and the yield for each t] c decay channel, the 16 
7T° recoil mass distributions are fitted simultaneously with a binned maximum likelihood 
method. A Breit-Wigner function convolved with the instrumental resolution is used to 
describe the signal shape. An efficiency correction is not needed because of the small h c 
width and the good ir° mass resolution. The resolution function is channel- dependent and 
is obtained from MC simulation. The parameters M(h c ) and T(h c ) of the Breit-Wigner 
function are constrained to be the same for all 16 channels, which is essential for the decay 
modes with low statistics. For the recoil mass fit to each channel, the background shape 
is obtained from the r\ c mass side bands (2300—2700, 3070— 3200 MeV/c 2 ), and the signal 
and the background normalizations for each mode are allowed to float. The summed and 
mode-by-mode fit results are shown in Figs. [1] and |2j respectively. The \ 2 P er degree of 
freedom for this fit is 1.60, where sparsely populated bins are combined so that there are 
at least seven counts per bin in the x 2 calculation. The parameters of the h c resonance are 
determined to be M(h c ) = 3525.31 ± 0.11MeV/c 2 and F(h c ) = 0.70 ± 0.28 MeV, where the 
errors are statistical only. 

The MC-determined selection efficiency e« and yield iVj for each r] c decay mode are 
listed in Table [Til Based on these numbers, we can calculate the product branching ra- 
tios i3i(^(3686) — > 7r°h c ) x B<i{h c — > 7?7 C ) x B^,{rj c — > Xi). The branching ratio for r] c — > X-i 
for each of the 16 final states Xi can then be obtained by combining our measurements 
with Si 0(3686) ->• u°h c ) x B 2 (h c ->- jt] c ) = (4.36 ± 0.42) x 10" 4 , the average of two recent 
measurements by CLEO [ll[ and BESIII [13]. These branching ratios, with both statistical 
and systematic errors, are presented in Section IVT1 

B. Measurement of r\ c resonance parameters 

In addition to determining the h c resonance parameters, we can also measure the r] c 
mass and width with the same event sample. The decay chain h c — > 77/c r] c — > Xi is 
reconstructed and kinematically fitted in the 16 r] c final states Xi. For candidates with 
satisfactory kinematic fits, we use the resulting track and photon momenta to compute the 
hadronic mass. We populate distributions of this hadronic mass by removing our previous 
El photon-energy and M(t] c ) requirements and selecting candidates inside a 7r° recoil mass 
window of ±5MeV/c 2 around the h c mass, keeping all other criteria unchanged. 

The line shape for the r\ c signal for these fits is parameterized as (E^ x BW(m) x fd(Ey))<g) 
Ri(m), where BW(m) is the Breit-Wigner function for rj c as a function of the invariant mass 

m of the decay products for each channel, E^(m) = M ^ ft ~) m is the energy of the transition 

photon in the rest frame of h c , and fd^E^) is a function that damps the divergent tail due to 
the E^ factor, which incorporates the energy dependence of the El matrix element and the 
phase-space factor. Ri(m) is the signal resolution function for the ith decay mode, which 
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7i° recoil mass (GeV/c 2 ) 

FIG. 2: The simultaneously fitted ir° recoil mass spectra in -0(3686) — > TT°h c ,h c —¥ "yr] c , r\ c — > Xi 
for the 16 final states Xi. 



is parameterized by double Gaussians to account for the distorting effects of the kinematic 
fit and detector smearing. The damping function that we use was introduced by the KEDR 
collaboration |3C 



fd{E 1 



El 



E y E + (E 1 — E ) 



where Eq = E^{m Vc ) is the i?l-transition-photon peak energy. The ?7 c -candidate hadronic 
invariant mass spectra from low and high side bands in the h c mass (3500—3515, 
3535— 3550 MeV/c 2 ) are used to obtain the background functions for the r\ c mass fit. To 
mitigate the effects of bin-to-bin fluctuations, these side-band mass spectra are smoothed 
before fitting. A toy MC study was performed to test the effect of the smoothing and it 
was demonstrated to be a robust procedure that does not systematically distort the fit re- 
sults. The channel-by-channel signal and background normalizations are free parameters 
determined by the fit. 

We ignore the effect of interference between the signal and background, which was con- 
sidered in the previous measurement of -?/>(3686) — > 7% 20|, because the branching ratio of 
h c — > IVc is about 50% (branching ratio of Ml transition ^(3686) — > •yr] c is about 0.3%). 
The radiative decay of h c — > 70~ should be the same level of -?/>(3686) — > 70~, in this case, 
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TABLE II: MC-determined efficiencies e, and yields Ni for ^(3686) — > n°h c ,h c — > ^r] c ,r] c Xi, 
where Xi refers to the 16 final states . 



iviocie 




AT. 


PP 


zz.z 


1 rr o 1 >| r 
10. o ± 4.0 


TT^TT TT^TT 


in /; 

Iz.o 


100. ± 11. 6 


I\- IX xv iv 


D.D 


(!(! J. OK 
D.D it Z.O 


Pl J\ TT TT 


o. f 




PpTT+TT 


7 a 
f .0 


1 Q fl _l_ K A 

iy.u ± 0.4 


TT+1T 1T+1T TT TT 


0.4 


ou.o ± y.u 


K+ K TT+TT TT IT 


7 
Z. f 


1 n Q J- A o 

lU.o ± 4.y 


K + K~TT 


114 


U Q J- Q O 

o^t.y in y.z 


ppIT 


8.9 


14.4 ± 4.6 


K^K+TT* 


8.9 


107.1 ± 11.8 


K^K+tt+tt+tt^ 


3.4 


43.3 ± 8.0 


TT + lT~ri 


4.3 


32.9 ± 6.7 


K+K-r] 


3.0 


6.7 ± 3.2 


ir + ir^ ir + ir~ r] 


1.9 


38.6 ± 7.6 


7r + 7r _ 7r 7r 


5.5 


118.4 ± 12.8 


7T + 7T~ 7T + 7T~ 7T 7r° 


2.2 


175.2 ± 17.3 


Total 




831.9±35.0 



the non-?7 c intensity in h c is much smaller than that for f/>(3686) — > , ~fri c . 

Figs. [3] and H] show the hadronic-mass-fit results. The r] c mass and width are determined 
to be M(r] c ) = 2984.49 ± 1.16MeV/c 2 and T(r] c ) = 36.4 ± 3.2 MeV, where the errors are 
statistical. The x 2 P er degree of freedom for this fit is 1.52, using the same x 2 calculation 
method to accommodate low-statistics bins as for the fit to the tt° recoil mass spectrum. 

V. SYSTEMATIC UNCERTAINTIES 



A. h c parameter measurements 

The systematic uncertainties for the M(h c ) and T(h c ) measurements are summarized 
in Table III II All sources are treated as uncorrelated, so the total systematic uncertainty is 
obtained by summing them in quadrature. The following subsections describe the procedures 
and assumptions that led to these estimates of the uncertainties. 

1. Energy calibration 

The potential inconsistency of the photon-energy measurement between data and MC is 
evaluated by studying -^(3686) — > jXci,2 ( Xci,2 — > 7^/0 > J/^ ~ > ^"V") f° r photons with low 
energy and radiative Bhabha events for photons with high energy. Discrepancies of 0.4% 
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M(hadrons) (GeV/c 2 ) 




M(hadrons) (GeV/c 2 ) 



FIG. 3: (a) The hadronic mass spectrum in -0(3686) — > ir°h c , h c — > r fr] c , rj c — > Xi summed over the 
16 final states JQ. The dots with error bars represent the hadronic mass spectrum in data. The 
solid line shows the total fit function and the dashed line is the background component of the fit. 
(b) The background-subtracted hadronic mass spectrum with the signal shape overlaid. 
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M(hadrons) (GeV/c 2 ) 

FIG. 4: The simultaneously fitted hadronic mass spectra for the 16 r\ c decay channels. 



in the energy scale and 4% in the energy resolution between data and MC are found. We 
vary the photon response accordingly and take the changes in the results as the estimated 
systematic error. For the M(h c ) measurement, besides the above studies, the reconstructed 
photon position and error matrix are taken into account as additional sources of uncertainty. 



2. Signal shape 

The uncertainty associated with the h c signal shape in the tt° recoil mass spectrum 
includes contributions from the photon line shape and the 1-C kinematic fit. We estimate 
these by determining the changes in results after reasonable adjustments in the photon 
response. The photon-energy resolution is estimated with the control sample ^(3686) — > 
7X c2 . As above, the energy resolution in data is found to be about 4% worse than in the 
MC simulation. We correct for this discrepancy by adding single- Gaussian smearing to the 
energy of the 7r° daughter photons and then using the alternative 7r° shape to redo the fit. 
The changes in results are assigned as the systematic errors. 
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TABLE III: The systematic errors for the h c mass and width measurements. 



Sources 


AM hc (MeV/c 2 ) 


AT hc (MeV) 


Energy calibration 


0.13 


0.07 


Signal shape 


0.00 


0.06 


r nnng range 


n da 


U. ID 


Binning 


0.02 


0.01 


Background shape 


0.01 


0.08 


Background veto 


0.01 


0.08 


Kinematic fit 


0.03 


0.03 


Mass of ^(3686) 


0.03 


0.02 


Total 


0.14 


0.22 



3. Fitting range and binning 

The systematic uncertainties due to the fitting of the ir° recoil mass spectrum are evalu- 
ated by varying the fitting range and the bin size in the fit. The spreads of results obtained 
with the alternative assumptions are used to assign the systematic errors. 



4- Background shape 



To estimate the uncertainty associated with the side-band method for assigning back- 
ground function shapes, we use an ARGUS function [3l| as an alternative background de- 



scription for each channel and record the changes in the fit results. 



5. Background veto 



The systematic uncertainties associated with the requirements to suppress background 
are estimated by varying the excluded ranges. 



6. Kinematic fit 



Systematic uncertainties caused by the kinematic fit are studied by tuning the tracking 
parameters and error matrices of charged tracks and photons based on the data. Control 
samples of J/vp -» (ft / n (980), <ft ->■ K + K~, / (980) -)■ 7r + 7r~, and ?/>(3686) -)■ 7x c j are used 
for this purpose j32[. Channel-by-channel changes of M(r) c ) and T(r] c ) are calculated after 
the tuning and then averaged by yields and taken as systematic errors. 



7. ?ft(3686) mass 

The systematic uncertainties of the M(h c ) and T(h c ) determinations associated with the 
uncertainty in the ^(3686) mass are estimated to be 0.03 MeV/c 2 and 0.02 MeV, respectively. 
These are found by shifting M^( 3686 ) by one standard deviation according to the PDG value 
[2| and redetermining the results. 
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B. r/ c branching ratio measurements 



The systematic errors in the 7] c branching ratio measurements are listed in Tables IIV1 
All sources are treated as uncorrelated, so the total systematic uncertainty is obtained 
by summing them in quadrature. The following subsections describe the procedures and 
assumptions that led to the estimates of these uncertainties. 

TABLE IV: The systematic errors (in %) in the r\ c branching ratio measurements of the r\ c exclusive 
decay channels. 



Sources 


PP 2 


7T+7T 


) 2{K + K- 


) K + K 7T+7T 


ppir + n 


3(7T + 7T 


) K + K-2{tt+tt- 


) K+K-TT 


N(V>(3686)) 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


Tracking 


4.0 


8.0 


8.0 


8.0 


8.0 


12.0 


12.0 


4.0 


PID (K° s ) 


2.0 


6.0 


6.0 


4.0 


4.0 


8.0 


8.0 


2.0 


Photon eff 


3.0 


3.0 


3.0 


3.0 


3.0 


3.0 


3.0 


5.0 


Fit range 


2.2 


1.2 


2.6 


2.9 


1.5 


5.3 


3.3 


2.7 


Bkg shape 


10.3 


2.5 


4.7 


0.9 


0.3 


0.2 


3.5 


2.8 


Signal shape 


2.3 


2.3 


2.3 


2.3 


2.3 


2.3 


2.3 


2.3 


KmFit eff. 


7.0 


6.3 


7.0 


8.8 


10.8 


7.3 


4.2 


2.0 


Bkg veto 


5.9 


5.5 


1.1 


0.6 


3.1 


2.3 


5.2 


1.7 


Cross feed 


0.0 


2.5 


0.0 


0.0 


0.0 


0.0 


0.0 


1.4 


T] c decay models 


0.0 


2.1 


3.7 


0.6 


2.5 


0.0 


3.0 


4.6 


r; c line shape 


0.7 


0.8 


0.6 


0.9 


0.6 


0.6 


0.6 


0.7 


Sum 


15.7 


14.8 


14.9 


14.1 


15.7 


18.0 


17.8 


10.6 



Sources 


ppTT 






7T+7T Tj 


K+K-r] 2 


7T+7T )r; 7T" 




2(7r+7T-7r°) 


N(^(3686)) 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


Tracking 


4.0 


8.0 


12.0 


4.0 


4.0 


8.0 


4.0 


8.0 


PID (K° s ) 


2.0 


1.0 


1.0 


0.0 


2.0 


0.0 


0.0 


0.0 


Photon eff 


5.0 


3.0 


3.0 


5.0 


5.0 


5.0 


7.0 


7.0 


Fit range 


7.7 


2.1 


1.5 


0.6 


6.0 


1.8 


0.6 


2.0 


Bkg shape 


0.1 


4.7 


4.7 


0.1 


5.9 


0.8 


3.3 


1.6 


Signal shape 


2.3 


2.3 


2.3 


2.3 


2.3 


2.3 


2.3 


2.3 


KmFit eff. 


6.8 


6.8 


7.3 


2.0 


1.2 


6.7 


2.4 


2.4 


Bkg veto 


3.7 


0.7 


2.8 


11.8 


5.4 


14.7 


12.8 


5.5 


Cross feed 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1.3 


0.0 


ij c decay models 


5.8 


2.5 


5.2 


5.5 


8.1 


0.0 


0.1 


0.5 


rjc line shape 


0.6 


0.6 


0.8 


0.8 


0.7 


0.7 


0.7 


0.7 


Sum 


14.8 


13.2 


17.0 


15.4 


15.3 


19.4 


16.4 


13.3 



1. Tracking and photon detection 

The uncertainty in the tracking efficiency is 2% per track and the uncertainty due to 
photon detection is 1% per photon [33|. MC studies demonstrate that the trigger efficiency 
for signal events is almost 100%, so that the associated uncertainty in the results is negligible. 
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2. PID and Kg reconstruction 

The systematic uncertainties due to kaon and pion identifications are determined to be 
2% in Ref. Q. We choose J/ip K*°K° S ,K*° Kn to evaluate the efficiency of K° s 
reconstruction. The 1% difference between data and MC is assigned as the systematic error 
due to this source. 

3. Kinematic fitting 

The systematic errors associated with kinematic fitting are estimated by using the control 
samples of ?/>(3686) — > 7r 7r°J/'i/> with J/ip decay to hadronic final states, which have similar 
event topology as ^(3686) — > 7i°h c ,h c — > / -fi] c . The average efficiency difference between 
data and MC, with the same \ 2 requirements in the h c selection, is taken as the systematic 
uncertainty. 

4- Cross- feed 

To evaluate the effect of cross-feed among the 16 signal modes, we use samples of 50,000 
MC events per mode. We find that r] c — > 2(7t + it~),t] c — > K + K~tt and r] c — > 7t + tt~it tt are 
contaminated by r\ c — > KgK ± ir T with levels of 2.5%, 1.4%, and 1.3%, respectively. These 
numbers are assigned as the systematic errors associated with cross-feed. For other channels, 
this contamination is found to be negligible. 

5. r/ c decay models 

We use phase space to simulate r\ c decays in our analysis. To estimate the systematic un- 
certainty due to neglecting intermediate states in these decays, we extract invariant masses 
of T] c daughter particles from -?/>(3686) —> ■yr] c ,ri c — > Xj. We analyze MC samples generated 
according to these invariant masses. To illustrate, Fig. [5] shows the invariant-mass distribu- 
tion comparison between the data and MC for the decay mode r] c — > K^K^Tt^ . In addition, 
for channels with low statistics and well-understood intermediate states, MC samples with 
these intermediate states were generated according to the relative branching ratios given by 
PDG. The spreads of the efficiencies obtained from the phase-space and alternative MC are 
taken as the systematic errors. 

6. r/ c line shape 

Because of the r/ c mass window requirement in our event selection, the line shape of r\ c 
could be a source of systematic error in the measurement. We vary the input r\ c resonant 
parameters by one standard deviation to estimate the uncertainty due to this source. 
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FIG. 5: The dots show the mass spectra for ^(3686) — > jr] c ,i] c — > K^K^ir^ in data, and the solid 
lines are the corresponding mass spectra from the MC simulation. 

C. r] c parameter measurements 

Systematic errors for the M(r] c ) and T(rj c ) measurements are summarized in Table |Vl 
All sources are treated as uncorrelated, so the total systematic uncertainty is obtained by 
summing in quadrature. The following subsections describe the procedures and assumptions 
that led to the estimates of these uncertainties. 



TABLE V: The systematic errors for rj c parameter measurements. 



Sources 


M(r? c ) (MeV/c 2 ) 


rfa c ) (MeV) 


Background shape 


0.36 


1.45 


Fitting range 


0.03 


0.33 


Resolution description 


0.10 


0.02 


Mass-dependent efficiencies 


0.11 


0.27 


Mass-dependent resolutions 


0.00 


0.01 


Kinematic fitting 


0.33 


0.76 


Fitting method 


0.11 


0.40 


Sum 


0.52 


1.74 



1. Background shape 

Our standard background shape is the smoothed h c side-band shape. To estimate the 
systematic uncertainty due to the background procedure, we change the smoothing level and 
technique, and vary the h c side-band ranges. The largest changes in results among these 
alternatives are assigned as the systematic errors. 
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2. Fitting range 

The systematic uncertainties due to the fitting range are estimated by considering several 
alternatives to the standard fitting range of 2.3-3.2 GeV/c 2 , 2.4-3.2 GeV/c 2 , 2.5-3.2 GeV/c 2 , 
2.6-3.2 GeV/c 2 , and 2.3-3.15 GeV/c 2 . The systematic uncertainties are assigned to be the 
largest differences between the standard fit results and those from the alternative ranges. 

3. Resolution description 

In order to estimate the systematic uncertainties associated with the detector-resolution 
description, we use MC signal shapes obtained by setting the i] c width to zero as alternatives 
to double Gaussians. The changes in fit results between these two methods provide the 
systematic errors. 

4- Mass-dependent efficiency and resolution 

Since the i] c signal spreads over a sizable mass range, the uncertainties due to the use 
of mass-independent efficiencies and resolutions need to be estimated. Mass-dependent 
efficiencies and resolutions are determined from MC simulation and used as an alternative 
to the default assumption, and the resulting differences are taken to be the systematic errors. 

5. Kinematic fitting 

The method to evaluate the systematic errors due to the kinematic fitting procedure and 
momentum measurement is the same as that in the measurement of the h c parameters. 

6. Fitting method 

Because we use the smoothed side-band shape to describe the background, the potential 
for bias due to the smoothing technique must be considered. This was investigated with 
a toy MC study. We start with a signal sample for each of the 16 channels selected from 
our standard MC to have the same statistics as data. A corresponding background sample 
for each channel is constructed from the mass side bands in data. The hadronic-mass 
distributions for these samples are then treated with a variety of smoothing procedures and 
fitted. The ranges in the fit results are used to set the systematic errors from this source. 

VI. SUMMARY AND DISCUSSION 

In summary, we have studied the process -^(3686) — > ir°h c followed by h c — > jrjc with an 
exclusive-reconstruction technique. Using a sample of 106 million -^(3686) decays we have 
obtained new measurements of the mass and width of the h c and r\ c charmonium resonances, 
and of the branching ratios for 16 exclusive r\ c hadronic decay modes. 
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The total yield of events, measured by fitting the 7r° recoil mass spectrum, is 832 ± 35 
events, where the error is statistical only. With these events we measure the mass and width 
of the h c ' 

M(h c ) = 3525.31 ± 0.11 ± 0.14 MeV/c 2 , and 
T(h c ) = 0.70 ± 0.28 ± 0.22 MeV, 
where the first errors are statistical and the second are systematic. These results are con- 



sistent with the results of a previous inclusive measurement by BESIII [131 ]: 

M(h c ) = 3525.40 ±0.13 ±0.18 MeV/c 2 , and 

T(h c ) < 1.44 MeV (at 90% confidence level). 

The branching-ratio results i3i(^(3686) — > ir°h c ) x B 2 (h c — > 7r/ c ) x B 3 (r/ C — > Xi) and 
&i(jlc are given in Table |VI| quoted with the statistical and systematic errors of 

this measurement and, for £>3, an additional systematic error associated with the input 
branching-ratio product B\ (ij; (3686) — > ir°h c ) x B 2 {h c — > 7??c)- Most of our B^r/c — > Xi) 
branching-fraction results are consistent with PDG values [2| , and several branching fractions 
are measured for the first time. 



TABLE VI: Bi (-0(3686) -> ir°h c ) x B 2 {h c jr] c ) x B 3 (t] c X t ) and 3 (r/ c X,) with sys- 
tematic errors. The third errors in B3 measurement are systematic errors due to uncertainty of 
£i(V(3686) n h c ) x B 2 (h c -> 7%). 





6ixB 2 x B 3 (xl(T 6 ) 


B, (%) 


B 3 in PDG (%) 


pp 


0.65 ± 0.19 ± 0.10 


0.15±0.04±0.02±0.01 


0.141±0.017 


7r + 7r~7r + 7r~ 


7.51 ± 0.85 ± 1.11 


1.72±0.19±0.25±0.17 


0.86±0.13 


K+K-K+K- 


0.94 ± 0.37 ± 0.14 


0.22±0.08±0.03±0.02 


0.134±0.032 


K+K-TT+TT- 


4.16 ± 0.76 ± 0.59 


0.95±0.17±0.13±0.09 


0.61±0.12 


pp7T + 7T~ 


2.30 ± 0.65 ± 0.36 


0.53±0.15±0.08±0.05 


<1.2 (at 90% C.L.) 


TT + 7T~7T + TT~7T + 7T~ 


8.82 ± 1.57 ± 1.59 


2.02±0.36±0.36±0.19 


1.5±0.50 


K + K~TT + TT~TT~TT~ 


3.60 ± 1.71 ± 0.64 


0.83±0.39±0.15±0.08 


0.71±0.29 


K+K-iT 


4.54 ± 0.76 ± 0.48 


1.04±0.17±0.11±0.10 


1.2±0.1 


PplT 


1.53 ± 0.49 ± 0.23 


0.35±0.11±0.05±0.03 






11.35 ± 1.25 ± 1.50 


2.60±0.29±0.34±0.25 


2.4±0.2 


ifg J P£T ± 7r =F 7r ± 7r = F 


12.01 ± 2.22 ± 2.04 


2.75±0.51±0.47±0.27 




7T + 7r _ 7) 


7.22 ± 1.47 ± 1.11 


1.66±0.34±0.26±0.16 


4.9±1.8 


K+K-r) 


2.11 ± 1.01 ± 0.32 


0.48±0.23±0.07±0.05 


<1.5 (at 90% C.L.) 




19.17 ± 3.77 ± 3.72 


4.40±0.86±0.85±0.42 




7r + 7r _ 7r°7r° 


20.31 ± 2.20 ± 3.33 


4.66±0.50±0.76±0.45 




7T + 7r _ 7r + 7r _ 7r 7r° 


75.13 ± 7.42 ± 9.99 


17.23±1.70±2.29±1.66 





Combining our measurement of M(h c ) with the previously-determined mass of the cen- 
troid of the 3 Pj states leads to 

A M h f = (M(1 3 P)) - M(1 1 P 1 ) = -0.01 ± 0.11 (stat.) ± 0.15 (syst.) MeV/c 2 , (1) 

consistent with the lowest-order expectation that the IP hyperfine splitting is zero. 
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The line shape of r\ c was also studied from the El transition h c — > 777c , and the measured 
resonance parameters are: 



M{rj c ) = 2984.49 ± 1.16 ± 0.52 MeV/c 2 , and 
T(r] c ) = 36.4 ± 3.2 ± 1.7 MeV. 
These results are consistent with the recent BESIII results from -^(3686) — > 7?7 C [201 ] : 

M(r/ C ) = 2984.3 ± 0.6 ± 0.6 MeV/c 2 , and 
r(r? c ) = 32.0 ± 1.2 ± 1.0 MeV; 



and B-factory results from 77 — > rj c and B decays [17j, |18|. Because of the larger -?/>(3686) 
data sample that will be coming from BESIII and the advantage of negligible interference 
effects, we expect that h c — > ^ri c will provide the most reliable determinations of the i] c 
resonance parameters in the future. 
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